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ABSTRACT 

Proton  Dose  Assessment  to  the  Human  Eye  Using  Monte  Carlo  N-Particle 
Transport  Code  (MCNPX).  (August  2006) 

David  Bernhardt  Oertli,  B.S.,  Colorado  State  University 
Chair  of  Advisory  Committee:  Dr.  John  R.  Ford 

The  objective  of  this  project  was  to  develop  a  simple  MCNPX  model  of  the 
human  eye  to  approximate  dose  delivered  from  proton  therapy.  The  calculated  dose 
included  that  due  to  proton  interactions  and  secondary  interactions,  which  included 
multiple  coulombic  energy  scattering,  elastic  and  inelastic  scattering,  and  non-elastic 
nuclear  reactions  (i.e.,  the  production  of  secondary  particles).  After  benchmarking 
MCNPX  with  a  known  proton  simulation,  the  proton  therapy  beam  used  at  Laboratori 
Nazionali  del  Sud-INFN  was  modeled  for  simulation.  A  virtual  water  phantom  was  used 
and  energy  tallies  were  found  to  correspond  with  the  direct  measurements  from  the 
therapy  beam  in  Italy.  A  simple  eye  model  was  constructed  and  combined  with  the 
proton  beam  to  measure  dose  distributions.  Two  treatment  simulations  were  considered. 
The  first  simulation  was  a  typical  treatment  scenario — where  dose  was  maximized  to  a 
tumor  volume  and  minimized  elsewhere.  The  second  case  was  a  worst  case  scenario  to 
simulate  a  patient  gazing  directly  into  the  treatment  beam  during  therapy.  Dose 
distributions  for  the  typical  treatment  yielded  what  was  expected,  but  the  worst  case 
scenario  showed  the  bulk  of  dose  deposited  in  the  cornea  and  lens  region.  The  study 
concluded  that  MCNPX  is  a  capable  platform  for  patient  planning  but  laborious  for 
programming  multiple  simulation  configurations. 
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CHAPTER  I 

INTRODUCTION  AND  BACKGROUND 

MCNPX  Background 

The  Monte  Carlo  method  describes  a  broad  area  of  mathematics  in  which  physical 
processes  or  other  phenomena  may  be  simulated  using  statistical  methods  through 
random  numbers.  One  of  the  earliest  examples  of  using  random  or  chance  events  to 
ascertain  truths  hidden  from  common  perception  was  demonstrated  by  Dr.  de  Buffon 
(Zaidi  2003).  In  1777  he  illustrated,  by  randomly  tossing  needles  onto  a  table  where 
parallel  lines  were  drawn,  that  one  could  determine  experimentally  the  value  of  pi.  It 
was  about  1944  when  Ulam  and  von  Neumann  (Metropolis  and  Ulam  1949)  coined  the 
name  ‘Monte  Carlo’  during  the  World  War  II  Manhattan  Project.  It  was  adopted  from 
the  well-known  city  in  Monaco  famous  for  its  games  of  chance.  Ulam  and  von 
Neumann  led  us  into  the  modem  age  of  Monte  Carlo  techniques  and  their  application 
with  modem  computers  was  first  realized  in  1949,  in  simulations  of  neutron  transport 
(Zaidi  2003). 

Monte  Carlo  codes  all  follow  the  same  general  formulation.  A  model  is  created  which 
is  as  similar  as  possible  to  the  real  system  of  interest,  the  model  then  is  related  to 
potential  interactions  within  that  system  based  on  known  probabilities  of  occurrence,  by 
random  sampling  of  probability  density  functions  (pdfs).  As  more  interactions  occur 
within  the  simulation,  the  quality  of  the  reported  average  behavior  of  the  model 
improves,  and  therefore  statistical  uncertainty  decreases.  General  components  that  are 
found  in  all  Monte  Carlo  simulations  include: 


This  thesis  follows  the  style  of  Health  Physics. 
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•  Probability  density  functions:  the  physical  system  modeled  must  be  related  to 
reality  through  a  described  set  of  probability  density  functions. 

•  Random  number  generator:  for  uniformly  distributing  values  on  the  unit 
interval. 

•  Sampling  protocol:  a  prescribed  method  for  sampling  applicable  probability 
density  functions. 

•  Scoring:  the  outcome  from  the  sampling  must  be  tallied  for  quantities  of  interest. 

•  Error  estimation:  Statistical  error  (variance)  must  be  calculated  in  light  of  the 
number  of  trials  or  other  relevant  quantities. 

•  Optimization  techniques:  methods  for  reducing  the  variance  which  also  reduces 
required  computational  time. 

The  code  MNCPX  possesses  all  of  these  characteristics,  serving  as  a  general  purpose 
Monte  Carlo  radiation  transport  code  (Pelowitz  2005).  It  has  progressed  as  a  code  to  the 
point  where  it  tracks  nearly  all  particles  and  energies.  It  is  the  latest  generation  of 
transport  code  derived  from  60  years  of  development  at  Los  Alamos  National 
Laboratory. 

MCNPX  development  began  in  1994  as  program  extension  of  MCNP4B  and  LAHET 
2.8 — with  the  goals  of  modeling  all  particles  at  all  energies;  improving  physics  models; 
extending  of  neutron,  proton,  and  photonuclear  libraries  up  to  150  MeV;  and  improving 
variance-reduction  and  data-analysis  techniques  (Pelowitz  2005).  The  program  also 
incorporated  cross-section  measurements,  benchmark  experiments,  deterministic  code 
improvements,  and  better  transmutation  code  and  library  tools. 
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The  initial  release  of  MCNPX  occurred  in  October  of  1997,  and  an  extensive  beta  test 
team  was  formed  to  continually  improve  the  code.  Largely  this  team  was  driven  by 
the  feedback  of  MCNPX  users.  To  date,  approximately  1500  users  in 
approximately  300  institutions  around  the  world  have  had  the  opportunity  to  provide 
developer  feedback.  This  process  has  tested  MCNPX  in  a  myriad  of  applications  and 
conditions,  proving  it  a  useful  and  reliable  tool.  Applications  for  the  code  are  broad  and 
constantly  increasing.  Some  examples  include  the  following: 

•  Investigations  for  accelerator  isotope  production  and  destruction  programs, 
including  the  transmutation  of  nuclear  waste; 

•  Research  into  accelerator-driven  energy  sources; 

•  Investigations  of  cosmic-ray  radiation  backgrounds  and  shielding  for  high 
altitude  aircraft  and  spacecraft; 

•  Accelerator-based  imaging  technology  such  as  neutron  and  proton  radiography; 

•  Design  of  shielding  in  accelerator  facilities; 

•  Activation  of  accelerator  components  and  surrounding  groundwater  and  air; 

•  Investigation  of  fully  coupled  neutron  and  charged-particle  transport  for  lower 
energy  applications; 

•  Comparison  of  physics-based  and  table-based  data; 

•  Charged-particle  tracking  in  plasmas; 

•  Charged-particle  propulsion  concepts  for  spaceflight; 

•  Detection  technology  using  charged  particles  (i.e.,  abandoned  landmines); 


Nuclear  safeguards; 
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•  Radiation  protection  and  shielding;  and 

•  Medical  physics,  especially  proton  and  neutron  therapy. 

The  extensive  testing,  use,  and  constant  improvement  of  MCNPX  has  yielded  a 
reliable  and  useful  tool  for  radiation  transport  computations.  MCNPX  has  been  widely 
recognized  as  capable  of  producing  accurate  and  precise  data  for  use  in  radiation  patient 
treatment  planning  (Newhauser  2005).  Simulations  that  were  at  one  point  thought  to  be 
incalculable  have  been  resolved  in  reasonable  simulation  times  through  code 
improvements  and  increased  computing  power.  Specifically,  MCNPX  can  now  take  into 
account  3-dimensional  geometry,  coulombic  energy  loss,  energy  straggling,  multiple 
coulomb  scattering,  elastic  and  inelastic  scattering,  and  inelastic  nuclear  reactions  (i.e., 
the  production  of  secondary  particles).  The  advances  in  Monte  Carlo  algorithms  have 
led  to  the  potential  use  in  radiation  treatment  planning,  however,  until  recently,  attempts 
have  not  been  made  to  demonstrate  this  capability  (Newhauser  2005).  Treatment 
planning  with  charged  particles  is  an  especially  important  area  because  of  the  intense 
local  energy  deposition. 

Proton  vs.  Conventional  Radiotherapy 

Conventional  radiotherapy,  which  applies  photon  beams,  is  the  most  frequently 
utilized  modality  for  treating  localized  tumors.  Tumor  control  is  achieved  through 
energy  deposition  from  a  high  energy  beam  to  a  localized  group  of  tumor  cells.  The 
challenge  in  radiotherapy  has  always  been  to  maximize  dose  to  the  neoplastic  growth, 
while  minimizing  radiation  induced  morbidity  in  surrounding  healthy  tissues. 
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In  past  years,  photon  therapy  has  progressed  as  a  science,  especially  in  the  realm  of 
imaging  and  treatment  planning.  To  date  the  most  advanced  methodology  is 
intensity-modulated  radiation  therapy  (IMRT),  which  delivers  higher  doses  to  target 
cells  as  compared  to  dose  delivered  to  surrounding  normal  cells.  IMRT  limits  dose  to 
healthy  tissues  by  applying  numerous  radiation  fields  of  varying  intensities  from 
different  directions.  Fig.  1  illustrates  a  typical  photon  IMRT  assembly  and  patient 
planning. 


adapted  from  A  T.  Vfllavtcenao, 
Duke  University  Me  dical  Center 


Fig  1:  Typical  photon  IMRT  planning  and  machine  configuration  (Liu  2005) 


While  IMRT  reduces  dose  to  healthy  tissues,  it  requires  a  larger  treatment  volume — 
while  energy  deposition  to  the  tumor  area  increases,  it  requires  that  more  cells,  over  a 
greater  area  to  be  irradiated.  This  strategy  significantly  increases  the  integral  dose  to  the 
patient.  A  concern  for  this  type  of  treatment  is  the  potential  for  secondary  malignancies 
or  other  late  tissue  effects.  This  is  especially  a  concern  in  pediatric  patients — if  such 
patients  are  cured  of  their  primary  malignancies,  they  are  expected  to  have  relatively 
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long  lifespan  during  which  radiation-induced  effects  might  translate  into  health  risks 
later  in  life. 

In  comparison  to  photon  treatments,  protons  were  recognized  as  having 
potentially  superior  dose  distributions  as  early  as  1946  (Wilson  1946).  Protons,  like  all 
charged  particles,  exhibit  a  rapid  energy  loss  at  the  end  of  their  tracks.  The  resultant 
sharply  localized  peak  of  dose  is  known  as  the  Bragg  peak.  The  penetration  depth  of  the 
Bragg  peak  is  directly  related  to  the  initial  energy  of  protons  entering  a  target.  The 
greater  the  initial  energy,  the  deeper  the  occurrence  of  the  Bragg  peak.  A  comparison  of 
a  high  energy  photon  beam  and  a  modulated  energy  proton  beam  is  shown  in  Fig.  2. 


Depth 


Fig  2:  Dose/depth  comparison  of  photons  vs.  energy  modulated  protons  (Suit  2002) 

In  patient  treatment,  the  peak  energy  deposition  for  protons  is  achieved  by  selecting  a 
distribution  of  proton  energies,  where  their  combined  Bragg  peaks  result  in  a  uniform 
dose  over  the  target  volume.  Instead  of  exhibiting  a  single  Bragg  peak,  the  peaks  of 
various  proton  energies  combine  to  give  the  cumulative  effect  of  an  extended  Bragg 
peak,  maximizing  energy  distribution  across  the  patient  target  volume.  This  leads  to 
superior  dose  in  the  target  cells,  but  also  minimizes  exposure  to  healthy  tissues  beyond 


7 

the  target.  While  photons  will  continue  to  deposit  energy  beyond  the  target  volume, 
protons  deposit  all  their  energy  in  a  small  area — minimizing  exposure  beyond  the 
target  cells. 

In  some  cases  proton  therapy  has  overwhelming  advantages.  Ocular  tumors  are  an 
example.  Eye  function  depends  upon  a  number  of  tissues  in  close  proximity.  The 
cornea,  lens,  retina,  fovea,  and  optic  nerve  are  all  critical  organs  for  eyesight — 
radiotherapy  can  damage  these  organs  leaving  patients  without  eyesight.  However, 
proton  therapy  in  many  cases  has  shown  far  superior  results  compared  with  photon 
treatments  (Bertil  et  al.  2005).  Proton  treatment  is  often  the  first  choice  for  pediatric 
patients.  Since  cured  patients  have  potentially  long  life  spans,  radiation  exposures  are 
minimized  to  avoid  potential  radiation-induced  risks  later  in  life. 

Ocular  Tumors 

Uveal  melanoma  and  retinoblastoma  are  die  principal  tumors  originating  in  the  eyes 
of  adults  and  children,  respectively.  Melanoma  of  the  uveal  tract  is  the  most  common 
primary  intraocular  cancer  in  humans — accounting  for  about  12%  of  all  melanomas 
(Harbour  2003).  Malignant  melanomas  appear  most  commonly  in  the  choroid,  followed 
in  decreasing  order  by  the  ciliary  body,  the  iris,  the  conjunctiva,  and  the  skin  of  the 
eyelids.  The  rarest  of  the  melanomas  are  those  that  originate  from  the  cornea  or  the 
orbit.  Incidence  rates  for  these  melanomas  have  been  linked  to  genetic  predisposition, 
immunological  or  hormonal  alteration,  and  environmental  and  occupational  factors,  such 
as  indoor  working  conditions,  exposure  to  chemicals,  radiofrequency  radiation  and 
ultraviolet  light. 
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The  appearance,  size  and  location  of  uveal  melanomas  are  highly  variable,  but 
clinical  examination  and  ancillary  testing  can  accurately  diagnose  most  uveal  and 
iris  melanomas.  This  usually  occurs  through  the  use  of  ultrasonography  because 
melanomas  typically  have  lower  reflective  properties  than  surrounding  tissues. 
Fluorescein  angiography,  magnetic  resonance  imaging,  or  fine  needle  aspiration  biopsy 
are  also  possible  diagnostic  tools.  Various  tumors  are  illustrated  in  Fig.  3. 


Fig  3:  Uveal  melanomas  may  include  the  iris  (top  left),  choroid  (top  right) — and 
progress  outward,  or  they  may  extend  internally  (bottom  left).  (Finger,  Char,  Kimmel 

Cancer  Center  2006) 


Ocular  Tumor  Treatments 

To  date  the  most  common  treatment  for  uveal  melanomas  is  enucleation — removal  of 
the  entire  eyeball.  This  is  most  considered  mostly  in  cases  that  involve: 
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•  Large  tumor  size; 

•  Tumor  invasion  of  the  optic  nerve  head; 

•  Lack  of  access  for  other  treatment  options; 

•  Inability  to  return  for  follow-up;  and 

•  Patient  choice. 

In  the  1970s,  Zimmerman  and  his  colleagues  were  the  first  to  document  the 
phenomena  that  enucleation  actually  has  the  potential  to  speed  metastatic  progression 
(Zimmerman  et  al.  1978).  This  may  occur  because  additional  tumor  cells  are  released 
into  the  bloodstream  during  the  handling  of  the  eye  in  the  enucleation  process.  It  has 
also  been  observed  in  animals  that  primary  tumors  may  produce  angiostatin,  a  tumor 
suppressant  that  effectively  inhibits  angiogenesis  and  additional  metastatic  growth.  This 
study  which  has  been  repeated  many  times,  leads  some  to  the  conclusion  that  in-place 
treatment  is  preferable  when  possible. 

Perhaps  slightly  more  conservative  than  enucleation  is  local  resection  within  the  eye. 
When  a  patient  must  retain  vision  in  a  seeing  eye,  or  the  tumor  cells  are  in  the  only 
remaining  eye,  local  resection  might  be  a  desired  approach.  This  procedure  entails 
surrounding  the  area  of  resection  with  photocoagulation  barrages  over  several  weeks  to 
cut  off  the  choroidal  blood  supply  to  the  area.  The  portion  of  the  sclera,  choroid,  or 
retina  with  the  tumor  is  then  excised,  and  the  wall  defect  is  repaired  with  a  sclearal  graft. 
This  is  a  complicated  surgical  procedure,  which  has  a  high  probability  for  seeding  tumor 
cells  into  the  wound  or  outside  the  eyeball  (Albert  and  Polans  2003). 

Transpupillary  Thermotherapy  was  developed  as  a  less  evasive  method  for  treating 
intraocular  tumors.  The  most  current  treatments  are  applied  using  an  infrared  diode  laser 
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(projecting  a  spot  size  of  2-3  mm),  used  for  about  1  minute  per  exposure.  The  laser 
heats  lesions  to  a  temperature  great  enough  to  induce  necrosis  without 
photocoagulation  in  target  cells.  Limitations  of  this  treatment  include  depth  of 
treatment  at  3.9  mm,  potential  for  unfavorable  side-effects  such  as  macular  traction, 
retinal  vascular  occlusion,  macular  edema,  retinal  or  vitreous  hemorrhage,  and  visual 
field  defects.  Thermotherapy  is  usually  applied  concurrently  with  some  form  of 
radiotherapy  (Godfrey  1999). 

Radiotherapy  is  typically  applied  in  one  of  three  forms:  conventional  external-beam, 
brachytherapy,  and  charged  particles.  Conventional  external-beam  therapy  is  available, 
but  generally  not  applied  due  to  complications  discussed  previously.  Brachytherapy  is  a 
more  recent  development  aimed  at  delivering  localized  radiation  with  fewer  side  effects 
than  external  beam  therapy.  A  typical  configuration  involves  attaching  radioactive  seeds 
to  a  lead  or  gold  “plaque”  that  is  sewn  to  the  sclera  overlying  the  imbedded  tumor. 
Common  radionuclides  include  iodine-125,  ruthenium-106,  and  palladium- 103. 
Typically  80-100  Gy  is  delivered  to  the  tumor  apex  over  a  4-5  day  period,  with  a  2-mm 
margin  around  the  tumor  base  (Robertson  1983).  Limitations  include  the  invasive  nature 
of  treatment,  and  radiation  side-effects  including:  cataracts,  retinopathy,  papillopathy, 
and  neovascular  glaucoma  These  effects  are  strongly  dose-dependent  and  begin  to 
increase  sharply  at  40  Gy  (Parsons  1996),  which  is  below  the  minimum  level  needed  for 
local  tumor  control. 

As  discussed  previously,  charged  particle  therapy  is  a  means  for  delivering  highly 
focused  radiation  to  tumor  areas.  This  is  arguably  one  of  the  best  options  for  treatment. 
Zimmerman  and  his  co-workers  (Zimmerman  et  al.  1978)  identified  an  important  trend 
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in  metastatic  progression  following  invasive  treatments.  Enucleation  and  local  resection 
are  the  only  treatment  options  that  introduce  the  risk  factor  that  tumor  cells  may 
seed  in  other  tissues  as  a  result  of  the  invasive  handling  of  lesions.  Brachytherapy, 
or  the  use  of  radiotherapy  plaques,  is  also  fairly  invasive  and  introduces  surgical 
difficulties  that  affect  accuracy  of  plaque  placement.  Thermotherapy  is  less  invasive, 
but  not  as  effective.  The  limited  depth  of  treatment  is  a  significant  shortcoming,  and  it 
requires  concurrent  treatment  with  other  modalities.  Charged  particle  therapy,  or 
specifically  proton  therapy  does  have  limitations — it  can  be  cost  prohibitive,  it  has 
higher  potential  for  anterior  segment  complications,  and  neovascular  glaucoma 
occurrence  is  higher  than  in  other  treatments.  However,  proton  therapy  boasts  superior 
tumor  control  after  treatment  compared  with  other  radiotherapy  methods,  which  greatly 
lowers  the  risk  of  metastic  disease.  The  non-invasive  nature  of  treatment  is  also  a 
positive  aspect  of  proton  therapy.  Less  manipulation  of  the  eye  translates  into  less 
potential  for  unintended  spreading  of  tumor  cells,  and  accidental  changes  in  position  of 
the  treatment  target  volume.  In  spite  of  the  limitations,  proton  therapy  is  now  a  leading 
treatment  modality  for  ocular  tumors  (Levin  et  al  2005). 
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CHAPTER  II 
MAKING  OF  A  MODEL 


Benchmarking 

To  establish  that  the  MCNPX  code  was  functioning  as  expected,  a  known  simulation 
was  applied  locally.  A  10-cm  diameter  200-MeV  proton  beam  was  constructed  to 
irradiate  a  cylindrical  water  phantom.  Three  smaller  cylinders  were  placed  in  the  path  of 
the  beam  (inside  the  water  phantom)  to  test  the  qualitative  interactions  with  the  proton 
beam.  The  first  cylinder  was  defined  as  a  zero-importance  area  for  protons — meaning 
that  any  protons  that  enter  the  volume  ceased  to  exist.  The  second  cylinder  in  the  path  of 
the  beam  was  constructed  of  iron,  and  the  third  was  defined  as  a  void  space. 

Qualitatively  speaking,  one  would  expect  that  particles  would  cease  to  exist  in  the  zero- 
importance  area,  and  that  a  “shadow”  would  be  observable  down-stream  of  the  source 
from  that  volume.  In  the  iron  cylinder  one  would  expect  greater  interaction 
corresponding  to  a  much  higher  proton  stopping-power  and  short  proton  path  length. 
Finally,  in  the  void  volume,  we  would  expect  to  see  no  interactions— and  a  greater  path 
length  for  particles  downstream  of  the  void  (due  to  a  decrease  in  particle  interactions  in 
their  path  of  travel  due  to  the  void). 

To  illustrate  the  results  of  this  benchmark,  a  mesh  tally  plot  was  used.  Mesh  tallies 
were  used  for  tracking  fluxes,  heating,  source  locations,  doses  and  other  tally  quantities. 
For  this  case,  a  heating  mesh  tally  was  applied  and  plotted.  The  results  are  shown  in  Fig. 


4. 


13 


Void 

Iron  Inclusion 

Zero-importance 

Region 


Fig  4:  Heating  mesh  tally  plot  for  a  200-MeV  proton  beam 

Results  of  the  2-D  plot  of  the  water  phantom  with  inclusions  yielded  what  one  would 
expect.  The  blue  color  scheme  indicated  cooler  regions,  while  the  red  showed  higher 
temperatures.  The  zero- importance  area  (far  left  inclusion)  indicated  desired  qualitative 
results.  Very  limited  interactions  were  found  inside  the  volume  (zero-importance  was 
applied  to  protons,  not  all  particles).  There  was  also  an  abrupt  change  in  proton 
interactions  downstream  of  the  zero-importance  area  as  evidenced  by  a  clear  shadow 
following  the  zero-importance  volume.  The  iron  inclusion  showed  greater  heating  as 
expected,  and  downstream  from  the  iron  inclusion,  the  length  of  travel  for  protons  was 
much  less  than  what  was  found  elsewhere.  The  last  inclusion,  a  void,  also  yielded  the 
expected  results.  There  were  no  interactions  in  the  void,  and  there  was  clearly  a  longer 
path  of  travel  for  protons  that  traveled  downstream  from  the  void.  A  flux  tally  was  also 
plotted  and  appropriate  values  were  obtained.  These  tests  confirmed  that  the  MCNPX 
code  was  simulating  the  proton  beam  at  least  qualitatively. 
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To  benchmark  a  beam  for  radiotherapy  simulation,  the  proton  therapy  beam  used  at 
Laboratori  Nazionali  del  Sud-INFN  was  selected  as  a  model  for  the  simulation. 

Built  in  collaboration  with  the  University  of  Catania  (Italy),  the  62-MeV  proton 
beam  is  used  for  treatment  of  ocular  pathologies  like  uveal  melanomas,  choroidal 
hemangiomas,  conjunctiva  melanoma,  eyelid  tumors,  and  embryonal  sarcomas.  Before 
each  treatment,  the  beam  is  calibrated  using  a  water  phantom  with  embedded  PTW 
34045  Markus  Ionization  Chambers.  The  ion-chambers  have  an  electrode  spacing  of 
1mm  and  sensitive  air-volume  of  0.02  cm3.  Dose  measurements  using  the  phantom  and 
embedded  ion-chambers  have  been  compiled  and  published  (Cirrone  2004). 

To  simulate  the  62-MeV  beam,  the  original  MCNPX  simulation  source  energy  and 
beam  diameter  were  adjusted  to  match  the  Laboratori  Nazionali  del  Sud-INFN  proton 
beam.  The  MCNPX  water  phantom  was  subdivided  into  progressive  layers  of  water  in 
which  deposited  dose  was  tallied.  This  enabled  measurement  of  energy  deposition  as  a 
function  of  depth  in  the  phantom.  The  simulation  yielded  the  expected  Bragg  peak 
behavior.  The  maximum  penetration  depth  is  slightly  greater  in  the  simulation.  This  is 
primarily  due  to  electron  straggling  effects,  and  lower  energy  cut-off  in  simulation.  The 
MCNPX  code  does  not  track  scattered  electrons  with  an  energy  value  lower  than  1  keV. 
This  approximation  eliminates  some  of  the  dose  due  to  delta-rays  in  the  target.  Likewise 
the  lower-energy  cut-off  for  protons  is  1  MeV,  which  lowers  the  overall  dose  in  the 
target  and  explains  the  slight  variation  between  the  two  sets  of  data.  Measured  values 
obtained  in  Italy  and  the  MCNPX  simulation  results  can  be  compared  in  Fig  5. 
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Fig  5  :  Measured  activity  in  water  phantom  (left),  MCNPX  simulation  results  (light) 


For  applying  treatment  to  individuals,  it  is  common  practice  to  introduce  a 
modulating  device  and  a  range  shifter  to  effectively  spread  out  the  Bragg  peak  and  vary 
the  maximum  distance  of  the  beam.  The  modulator  acts  on  the  beam  to  effectively  vary 
the  beam  energy  such  that  a  carefully  selected  spectrum  of  proton  energies  emerge  so 
that  a  uniform  dose  over  the  volume  of  interest  will  result.  This  is  usually  achieved  by 
placing  a  spinning  disk,  or  a  modulator  wheel,  between  the  source  and  target.  The 
spinning  wheel  has  different  material  characteristics  at  different  points  on  the  disk.  As 
the  proton  beam  interacts  with  the  modulator  material,  beam  characteristics  are  changed, 
creating  a  spectrum  of  proton  outputs.  Fig.  6  is  a  photo  of  a  modulator  wheel  used  in 
proton  therapy  at  Laboratori  Nazionali  del  Sud-INFN. 
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Fig  6:  Modulator  wheels  used  with  the  Laboratori  Nazionali  del  Sud-INFN  proton  beam 

(Cirrone  2005) 


If  one  sums  together  the  various  Bragg  peaks  of  individual  proton  energies  exiting  a 
modulation  wheel,  a  spread  out  Bragg  peak  can  be  obtained,  as  illustrated  in  Fig.  7. 


Fig  7:  A  spread  out  Bragg  peak  formed  from  the  sum  of  individual  Bragg  peaks 

(Kooy  2003) 


The  resultant  beam  deposits  much  higher  energies  over  a  larger  target  volume.  This 
allows  high-energy  treatment  to  a  tumor  volume  instead  of  a  tumor  point,  and  at  the 
same  time  minimizes  dose  beyond  the  tumor  cells. 


17 


To  model  a  modulator  wheel  in  MCNPX  is  a  challenging  task.  A  limitation  of  most 
Monte  Carlo  simulations  (MCNPX  included)  is  that  material  properties  cannot  vary 
with  time  during  a  simulation.  To  model  a  rotational  wheel  in  the  geometry  of  a 
simulation  is  impossible.  A  varying  source  however  is  possible,  but  to  find  the  energy 
spectrum  necessary  for  the  spread  out  Bragg  peak  is  a  not  trivial. 

The  MCNPX  code  is  able  to  vary  energy  emitted  from  a  source  if  an  energy  spectrum 
is  defined,  and  each  energy  level  is  weighted  as  to  how  often  it  is  emitted.  Trial  and 
error  was  the  first  method  for  finding  an  appropriate  spectrum,  but  this  was  quickly 
abandoned  because  of  the  complexities  involved.  To  obtain  precise  data  points,  a 
MATLAB  routine  was  made  that  incorporated  the  original  62-MeV  MCNPX  Bragg  peak 
data.  The  MATLAB  routine  combined  lower  energy  protons  into  the  spectrum  and  the 
frequency  with  which  each  should  be  introduced  into  the  beam  for  the  desired  effect. 

The  result  was  a  spread  out  Bragg  peak  that  closely  simulated  the  actual  data  collected 
from  the  Laboratori  Nazionali  del  Sud-INFN  proton  beam.  As  previously  discussed, 
slight  variations  between  the  two  can  be  attributed  to  lower  energy  cut-off  for  protons 
energies  less  than  1  MeV  and  electrons  less  than  1  keV.  A  comparison  of  the  two 
outputs  can  be  seen  in  Fig.  8. 
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Fig  8:  Measured  Spread  out  Bragg  peak  (left),  simulated  (right) 

With  this,  the  MCNPX  simulation  had  been  appropriately  benchmarked.  A  known 
simulation  was  used  to  find  expected  output — the  simulated  water  phantom  with  tally 
inclusions  calculated  results  as  expected.  After  selecting  the  Laboratori  Nazionali  del 
Sud-INFN  proton  beam  for  a  therapy  simulation,  the  previous  MCNPX  model 
characteristics  were  adjusted.  The  energy  level  and  beam  diameter  were  set  to  match 
conditions  used  with  the  Laboratori  Nazionali  del  Sud-INFN  beam.  Dose  distribution 
into  the  water  phantom  was  verified  with  actual  data,  and  a  spread  out  Bragg  peak  from 
the  simulation  matched  closely  the  actual  proton  beam  characteristics. 

The  Eye 

One  of  the  potential  benefits  of  proton  therapy  is  that  ocular  tumors  can  potentially  be 
treated  in-place,  avoiding  more  evasive  procedures  like  local  resection  or  enucleation. 
Preserving  eyesight  through  the  course  of  in-place  treatment  is  a  challenging  task 
because  the  critical  organs  necessary  for  eyesight  are  located  so  closely  together.  These 
critical  organs  for  eyesight  were  the  emphasis  of  the  MCNPX  model.  The  lens  of  the 
eye  has  long  been  understood  to  be  a  radiosensitive  organ.  Likewise,  radiation  exposure 
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to  the  cornea  can  cause  the  structure  to  become  opaque  and  obstruct  vision.  Since  many 
uveal  melanomas  appear  in  the  choroid  and  sclera  structures  of  the  eye,  these 
needed  to  be  modeled.  Lastly,  the  optic  nerve  bundle  near  the  posterior  axis  of  the 
eye  was  an  area  of  interest.  A  two  dimensional  rendering  of  the  eye  used  for  the 
MCNPX  modeling  is  shown  in  Fig.  9. 


Fig  9:  MCNPX  simulation  of  the  human  left  eye 


Dimensions  of  the  human  eye  were  adapted  from  Optics  of  the  Human  Eye  (Atchison 
and  Smith  2003).  The  model  was  constructed  using  concentric  spheres,  offset  as 
appropriate.  The  vitreous  humor  was  defined  to  be  the  volume  posterior  to  the  lens,  and 
the  anterior  chamber,  anterior  to  the  lens.  The  delineation  is  marked  in  Fig.  9  by  lines 
from  the  ends  of  the  lens  to  the  cornea.  The  optic  nerve  is  simulated  as  a  cylinder 
appropriately  offset  from  the  posterior  pole  of  the  eye.  The  subdivisions  in  the  outer 


wall  of  the  model  are  volumes  for  dosimetric  tallies. 
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Material  compositions  for  the  model  were  adapted  from  the  International 
Commission  on  Radiation  Units  and  Measurements  (ICRU)  Report  46,  Photon 
Electron,  Proton  and  Neutron  Interaction  Data  for  Body  Tissues.  This  report 
addressed  various  tissues  groups  in  the  body  and  defined  their  elemental  composition 
and  densities  for  purposes  of  radiation  dosimetry.  The  lens  of  the  eye  is  addressed 
directly  in  the  publication;  however  the  other  components  of  the  eye  were  approximated. 
Recent  studies  have  indicated  the  vitreous  and  anterior  humors  have  characteristics 
similar  to  the  properties  of  lymph  outlined  in  ICRU  46  (Macknight  2000  and  Brzezinski 
2004).  The  choroid  and  sclera  are  essentially  connective  soft  tissue  and  were  modeled 
correspondingly  as  soft  tissue.  The  optic  nerve  was  approximated  using  the  elemental 
composition  of  nerve  tissue  in  rats  (Stys  and  Lopachin  1996). 

Integrating  the  Eye  and  Beam 

With  the  eye  modeled  and  the  proton  beam  simulation  benchmarked,  the  next  step 
was  to  combine  the  two  together  in  a  realistic  way.  The  positioning  of  the  target  relative 
to  the  therapy  beam  is  of  course  vital  in  treatment  of  any  tumor.  For  treatment  of 
patients  with  proton  therapy,  this  is  typically  achieved  through  a  series  of  appointments 
at  the  treatment  facility.  Following  the  diagnosis  of  a  tumor,  the  size  and  dimensions  are 
identified  relative  to  other  parts  of  the  eye.  The  shape  and  location  of  the  tumor  are  then 
marked  by  an  ophthalmologist  using  tantalum  marker  clips  which  are  sewn  to  the  outside 
of  the  sclera  The  shape  and  position  of  the  target  volume  are  outlined  by  the  clips.  This 
form  of  marking  a  target  area  is  noninvasive  enough  that  the  clips  are  generally  not 
removed  from  the  patient  following  treatment. 
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One  to  two  weeks  after  the  insertion  of  surgical  clips,  the  patient  returns  to  the 
treatment  facility  to  establish  the  position  of  the  head  needed  for  the  proton 
treatment.  Two  main  tools  are  used  to  immobilize  the  patient’s  head,  a  custom 
block-bite  and  individual  mask.  The  block-bite  is  a  piece  of  material  the  patient  bites 
down  on  during  the  treatment.  It  is  custom  made  to  fit  the  individual  patient’s  teeth, 
ensuring  consistent  placement  in  the  patient’s  mouth  as  they  bite  down  on  the  block. 
Likewise,  the  individual  mask  is  a  form-fitting  mesh  that  fits  over  the  individual’s  face 
and  is  anchored  behind  the  head.  The  mask  is  set  firmly  against  the  patients  face  and 
head,  then  anchored  to  ensure  immobilization.  The  combination  of  the  block-bite  and 
custom  mask  allow  precise  and  accurate  placement  of  the  patient  relative  to  the  beam. 
One  treatment  facility  reports  reproducibility  of  1/10  mm  along  3  linear  axes  and  2 
rotational  axes  (Heufelder  2006).  Once  the  patient  has  been  effectively  immobilized  in  a 
reproducible  position,  a  series  of  diagnostic  x-rays  are  taken  to  precisely  measure  the 
positions  of  the  marker  clips  in  the  eye  relative  to  the  patients’  position  and  the  treatment 
beam. 

With  the  geometry  of  the  head  relative  to  the  beam  established,  the  treatment  team 
reviews  diagnostic  images  (MRIs,  CTs,  sonograms,  patient  placement  x-rays)  and  plans 
the  proton  target  dose  and  the  necessary  angle  of  the  eye  for  treatment.  To  fix  the  gaze 
of  the  patient  at  the  appropriate  angle,  a  small  light  or  other  focal  point  is  moved  into  a 
position  where  the  patient  can  focus  their  gaze  to  an  exact  point.  The  angle  of  the 
patient’s  gaze  can  help  to  minimize  the  dose  to  critical  organs.  For  example,  in 
treatment  of  a  posterior  melanoma,  it  may  be  possible  for  the  patient  to  fix  their  gaze 
looking  toward  the  ceiling,  effectively  moving  the  cornea  and  lens  out  of  the  proton 
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treatment  path.  After  combining  all  the  information  together  with  a  therapy  plan,  a 
feasibility  check  is  done  and  the  patient  is  ready  for  treatment.  A  typical  treatment 
configuration  is  shown  in  Fig.  10. 


Fig  10:  Typical  proton  treatment  configuration  (Heufelder  2006) 

To  combine  the  eye  and  the  beam  to  a  similar  configuration  for  the  MCNPX 
simulation,  two  cases  were  considered.  The  first  was  a  typical  treatment  scenario,  with 
the  intent  was  to  irradiate  a  tumor  area  at  the  posterior  of  the  eye  with  minimal  effect  to 
critical  organs  of  the  eye.  The  second  was  a  worst  case  scenario,  where  the  patient’s 
gaze  changes  to  place  the  cornea  and  lens  directly  into  the  path  of  the  beam. 
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CHAPTER  ID 

RESULTS  AND  DISCUSSION 

Two  cases  were  taken  for  proton  radiotherapy  simulation.  The  first,  a  typical 
scenario  where  exposure  to  the  eye  was  minimized  and  a  potential  tumor  volume  was 
exposed  to  maximum  proton  dose.  The  cancer  simulation  was  a  uveal  melanoma  in 
either  the  choroid  or  sclera  of  the  eye.  The  second  case  was  a  worst  case  scenario  in 
which  the  patient  changed  their  gaze  to  look  directly  into  the  proton  beam — maximizing 
dose  to  the  cornea  and  lens. 

Just  as  with  a  patient  proton  beam  treatment  configuration,  beam  adjustments  in  the 
simulation  needed  to  be  made  to  contour  the  beam  to  the  treatment  volume.  Two 
modifications  were  necessary  to  contour  the  beam  appropriately.  A  range  modifier  was 
introduced  to  adjust  the  maximum  penetration  depth  of  protons  from  the  beam,  and  a 
modulator  to  change  the  energy  spectrum  creating  a  spread  of  the  Bragg-peak  over  the 
appropriate  target  volume.  A  range  modifier  of  water  was  modeled  between  the  proton 
source  and  the  eye,  lowering  the  maximum  distance  traveled  by  individual  protons,  and 
limiting  delivered  dose  outside  the  target  volume.  Likewise  a  modulator  wheel  was 
simulated  to  spread  out  the  Bragg-peak,  effectively  distributing  proton  dose  over  the 
target  volume.  The  same  spread  out  Bragg  peak  simulated  in  the  proton  beam 
benchmark  was  used  for  the  therapy  simulations. 

At  first,  the  eye  model  was  used  with  all  material  properties  set  to  water.  This 
allowed  output  data  to  be  compared  qualitatively  with  the  benchmark  data  in  the  water 
phantom.  Results  yielded  dose  distributions  and  proton  penetration  depths  that  were 


consistent  with  the  benchmarked  results. 
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Next  the  materials  of  the  eye  were  changed  to  material  compositions  that  more 
closely  matched  the  actual  composition  of  the  eye  and  a  calculation  was  made  for 
typical  patient  treatment  scenario.  For  this  typical  treatment  scenario,  the  proton 
beam  was  positioned  to  give  minimal  dose  to  the  lens  and  cornea.  The  majority  of  the 
proton  dose  was  deposited  in  the  choroid/sclera  and  vital  organs  for  eyesight  were 
largely  spared. 

To  ensure  the  precision  of  these  results,  the  MCNPX  code  includes  ten  standard 
statistical  indices.  All  of  the  proton  beam  simulations  earned  passing  marks.  Two  of  the 
indices  are  discussed  here.  First,  the  tally  mean  measured  fluctuation  of  tally  counts  in 
relation  to  fluctuation  in  the  number  of  particles  simulated.  There  should  be  only 
random  fluctuations  in  tally  values  with  increased  particle  histories,  which  was  the  case 
with  the  proton  beam  simulations.  Another  statistical  measure  is  the  relative  error  (R). 
This  is  perhaps  the  most  scrutinized  of  the  statistical  indices.  It  relates  the  tally  mean 
with  the  overall  uncertainty — specifically  it  is  the  ratio  of  the  standard  deviation  of  the 
tally  mean  to  the  overall  mean.  R  is  typically  interpreted  against  the  criteria  shown  in 
Table  1  (Shultis  and  Faw  2006). 


Table  1.  Interpretation  of  the  relative  error  R. _ 

Simulation  R  Value _ Quality  of  Talley _ 

>  0.5  Meaningless 

0.2  to  0.5  Factor  of  a  few 

<  0.1  Reliable  (except  for  point/ring  detectors) 

<  0.05  Reliable  even  for  point/ring  detectors 


The  relative  error  for  the  typical  treatment  scenario  and  the  worst  case  scenario  were 
0.008  and  0.0012,  respectively — reflecting  extremely  low  relative  error.  After  numerous 
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simulations  it  was  clear  that  following  300,000  histories  was  sufficient  to  obtain 
excellent  statistics.  As  an  added  precaution,  these  calculations  were  made  several 
times  to  ensure  similar  results.  In  addition  to  tally  mean  and  relative  error,  300,000 
particles  simulations  attained  excellent  statistics  in  the  remaining  eight  statistical  indices 
in  MCNPX. 

To  obtain  a  graphical  representation  of  what  was  occurring,  MCNP  Visual  Editor 
Version  12N  was  used.  The  MCNPX  input  file  was  loaded  in  Visual  Editor  to  illustrate 
the  geometry  of  the  problem.  In  Fig.  11,  the  treatment  beam  travels  through  a  water 
range  modifier  and  is  incident  on  the  outer  structure  of  the  eye.  Arrows  have  been  added 
to  indicate  the  path  traveled  by  the  proton  beam. 


Fig  11:  Visual  Editor  illustration  of  the  typical  treatment  geometry 
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To  obtain  the  delivered  dose  to  various  parts  of  the  eye,  the  model  was  subdivided 
into  dosimetric  volumes,  which  include  the  lens,  cornea,  anterior  and  vitreous 
humors,  and  a  series  of  volumes  in  the  wall  of  the  eye  (see  Appendix  A).  These 
volumes  are  illustrated  in  Fig.  12. 


Fig  12:  Dosimetric  volumes  in  the  eye  model 


For  the  typical  treatment  scenario  the  volume  7R  was  identified  to  contain  the 
cancerous  tumor.  The  simulation  was  run  to  maximize  dose  to  this  portion  of  the  eye, 
and  minimize  it  elsewhere.  The  results  for  the  typical  treatment  scenario  are  found  in 


Table  2. 
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Table  2.  Dose  distribution  for  typical  treatment  scenario 


Dose  Volume 

Dose  per  Fraction  (Gy) 

Total  Dose  (Gy) 

Cornea 

0.60 

2.41 

Anterior  humor 

0.09 

0.37 

Lens 

0.09 

0.36 

Vitreous  humor 

4.93 

19.73 

Optic  Nerve 

1.06 

4.26 

1R 

11.01 

44.02 

2R 

10.86 

43.45 

3R 

10.25 

40.99 

4R 

9.71 

38.85 

5R 

9.67 

38.70 

6R 

10.51 

42.02 

7R 

12.50 

50.01 

8R 

10.22 

40.89 

9R 

1.56 

6.23 

10R 

0.02 

0.09 

1L 

0.02 

0.06 

2L 

0.01 

0.04 

3L 

0.01 

0.03 

4L 

0.00 

0.01 

5R 

0.00 

0.00 

6L 

0.00 

0.00 

7L 

0.00 

0.00 

8L 

0.00 

0.00 

9L 

0.00 

0.00 

10L 

0.00 

0.00 

A  typical  proton  therapeutic  dose  for  uveal  melanoma  in  this  portion  of  the  eye  is 
about  50  Gy  spread  over  four  fractions,  which  translates  into  four  treatments  of  12.5  Gy 
delivered  to  the  patient  (Heufelder  2006).  To  match  this  treatment  as  seen  in  Table  2, 
the  simulation  was  optimized  to  maximize  dose  to  7R  (12.5  Gy  per  fraction,  50.01  Gy 
overall)  while  minimizing  dose  elsewhere.  The  left  side  of  the  eye  received  doses  at 
least  three  orders  of  magnitude  lower  than  those  of  the  right  side.  The  optic  nerve  was  a 
concern,  not  because  the  tissue  is  as  radiosensitive,  but  because  it  was  directly  in  the 
path  of  the  proton  beam.  Results  showed  however,  that  dose  to  the  optic  nerve  was 
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relatively  small — only  about  1  Gy  per  fraction.  This  is  within  acceptable  the  limit  of  10 
Gy  to  the  optic  nerve  for  treatment  with  protons  (Jones  and  Errington  2000).  For 
each  treatment  the  cornea  received  less  than  3  Gy,  well  within  the  acceptable  limit 
of  15  Gy  (Simonva  2002).  For  the  lens  of  the  eye,  special  effort  is  made  in  therapy  to 
keep  doses  at  an  acceptable  limit — typically  less  than  8  Gy  (Jones  and  Errington  2000). 
In  this  simulation,  cumulative  dose  to  the  lens  of  the  eye  from  the  4  fractions  was  only 
0.36  Gy,  while  the  dose  to  the  target  volume  was  over  50  Gy.  The  dose  distribution  is 
shown  graphically  in  Fig  13. 


Fig  13:  Dose  distribution  for  typical  treatment  case 

For  the  worst  case  scenario,  the  beam  was  situated  to  treat  the  same  cancerous  target 
volume,  with  the  lens  of  the  eye  in  the  path  of  the  beam.  Since  the  object  of  this 
simulation  was  to  mimic  a  patient  gazing  into  the  beam  during  treatment,  the  same  dose 
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profile  from  the  typical  treatment  scenario  was  used.  The  geometry  for  this  simulation 
was  illustrated  using  Visual  Editor  and  arrows  were  drawn  in  to  represent  the  proton 
beam.  A  water  range  modifier  can  also  be  seen  at  the  bottom  of  Fig.  14. 


Fig  14:  Visual  Editor  illustration  of  the  worst  case  geometry 
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The  results  from  the  worst  case  simulation  are  shown  in  Table  3. 


Table  3.  Dose  distribution  for  worst  case  scenario _ 

Dose  Volume  Dose  per  Fraction  (Gy) _ Total  Dose  (Gy) 


Cornea 

9.02 

36.08 

Anterior  humor 

25.38 

101.52 

Lens 

38.72 

154.89 

Vitreous  humor 

5.96 

23.84 

Optic  Nerve 

0.00 

0.00 

1 R 

0.16 

0.64 

2R 

0.14 

0.57 

3R 

0.09 

0.37 

4R 

0.01 

0.05 

5R 

0.00 

0.01 

6R 

0.00 

0.00 

7R 

0.00 

0.00 

8R 

0.00 

0.00 

9R 

0.00 

0.00 

10R 

0.00 

0.00 

1L 

0.11 

0.42 

2L 

0.09 

0.35 

3L 

0.07 

0.28 

4L 

0.06 

0.25 

5R 

0.05 

0.21 

6L 

0.03 

0.14 

7L 

0.00 

0.02 

8L 

0.00 

0.00 

9L 

0.00 

0.00 

10L 

0.00 

0.00 

Surprisingly,  the  dose  to  the  cancerous  volume  in  the  eye  is  zero  for  this 
configuration.  The  majority  of  the  proton  energy  is  deposited  within  the  cornea,  anterior 
chamber  and  lens  of  the  eye.  If  one  fraction  of  therapy  was  conducted  in  this  manner, 
just  over  9  Gy  would  be  delivered  to  the  cornea  and,  if  this  configuration  occurred  for 
the  duration  of  treatment,  the  patient  would  experience  over  36  Gy  to  the  cornea  One 
would  expect  a  detrimental  effect  to  the  cornea  at  a  threshold  of  about  15  Gy  (Simonova 
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2002).  With  such  a  high  dose,  one  might  observe  keratitis  and  corneal  ulcerations 
possibly  leading  to  perforation  of  the  cornea.  Likewise,  the  lens  incurred  high 
doses.  For  a  single  fraction  in  this  position,  38.7  Gy  would  be  delivered  to  the  lens. 

For  the  treatment  series,  the  total  lens  dose  would  be  155  Gy,  compared  to  the  accepted 
tolerable  dose  of  8  Gy  (Jones  and  Errington  2000).  One  would  expect  severe  visual  loss 
due  to  the  lens  becoming  opaque  in  this  treatment  scenario.  The  cumulative  dose  in  the 
tumor  volume  for  this  configur  ation  is  zero,  which  is  illustrated  in  Fig  15. 


Fig  15:  Dose  distribution  for  worst  case  scenario 


The  data  from  the  worst  case  may  be  the  most  useful  of  the  simulation  data. 
Published  data  are  virtually  nonexistent  for  a  worst  case  scenario  like  this.  Of  course 
there  are  minimal  human  in  vivo  data  for  dose  reconstruction  of  this  type — this 
simulation  gives  an  approximation  that  could  be  used  to  reconstruct  dose  following  an 
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accident  in  treatment.  If  a  patient  was  to  alter  their  gaze  directly  into  the  treatment 
beam,  one  could  reconstruct  the  dose  by  the  fractions  of  time  in  each  treatment 
position.  The  fractions  of  dose  from  each  position  would  translate  to  dose  delivered 
to  the  lens  and  cornea,  and  dose  lost  from  the  target  volume. 
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CHAPTER  IV 
CONCLUSIONS 


Summary 

The  primary  objective  of  this  project  was  to  develop  a  model  of  the  human  eye  using 
the  computer  code  MCNPX  that  approximates  dose  delivered  during  proton  therapy. 

The  calculations  considered  proton  interactions  and  secondary  interactions,  which 
included  multiple-coulombic-energy  scattering,  elastic  and  inelastic  scattering,  and  non¬ 
elastic  nuclear  reactions  (i.e.,  the  production  of  secondary  particles).  After 
benchmarking  the  MCNPX  code  with  a  known  proton  simulation,  the  proton  therapy 
beam  at  Laboratori  Nazionali  del  Sud-INFN  was  modeled  for  simulation.  A  virtual 
water  phantom  was  used  and  energy  tallies  corresponded  with  the  direct  measurements 
from  the  therapy  beam  in  Italy.  An  eye  model  was  created  and  combined  with  the 
simulated  proton  beam.  Dose  was  tallied  for  two  specific  treatment  cases;  first,  a  typical 
treatment  case  where  dose  to  the  tumor  was  maximized,  while  doses  to  critical  organs 
for  eyesight  were  minimized.  The  second  case  was  a  simulation  where  the  patient  gazed 
directly  into  the  beam  during  treatment.  Doses  for  both  cases  were  calculated 
throughout  the  eye:  the  lens,  cornea,  anterior  and  vitreous  humors,  optic  nerve,  and  a 
series  of  volumes  in  the  wall  of  the  eye. 

Results  for  a  typical  treatment  scenario  yielded  the  desired  outcome.  The  dose  to  the 
tumor  volume  was  at  therapeutic  levels  and,  at  the  same  time,  doses  to  the  cornea,  lens 
and  optic  nerve  were  well  within  acceptable  limits.  For  the  worst  case  scenario,  the 
result  was  a  large  dose  to  the  cornea,  lens,  and  aqueous  humors.  The  tumor  did  not 
receive  any  dose  in  this  configuration. 
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MCNPX  Conclusions 

The  premise  for  this  endeavor  was  that  the  MCNPX  code  is  a  capable  platform, 
largely  under-utilized  for  this  type  of  application.  While  true,  there  maybe  some 
reasons  for  its  scarce  use  in  the  past,  e.g.,  the  complexity  of  data  input  and  limited 
flexibility  in  geometry  and  treatment  modifications.  The  user  interface  for  the  MCNPX 
code  is  simply  line  by  line  coding — an  interface  updated  little  with  the  progression  of 
technology  in  past  years.  Other  programming  languages  have  progressed  to  more  user 
friendly  interfaces;  like  Visual  Basic  or  Visual  C++,  characterized  by  mouse  controlled 
graphics,  drag-and-drop  icons,  and  common  code  progressions  that  can  be  duplicated 
and  inserted  into  new  programs.  Input  into  the  MCNPX  code  in  comparison  is 
laborious  and  antiquated.  During  the  construction  of  the  eye  model,  changes  in 
geometry  were  made  that  propagated  errors  throughout  the  rest  of  the  code.  A  single 
change  in  geometry  necessitated  careful  inspection  of  cell  definitions,  tally  boundaries 
and  volumes.  One  change  might  require  careful  inspection  of  three  or  four  sections  of 
the  code  to  ensure  all  variables  were  still  correctly  defined  and  compatible.  If  errors 
were  still  present  in  the  edited  code,  the  error  messages  generated  were  often  unable  to 
identify  the  specific  error,  but  instead  indicated  the  need  for  line  by  line  analysis  to 
restore  the  simulation  function.  Worse  yet,  there  were  a  handful  of  potential  errors  that 
did  not  generate  error  messages.  In  some  cases,  the  simulation  was  conducted  with 
obvious  errors  and  produced  erroneous  output.  Many  of  the  inconveniences  that  make 
the  MCNPX  code  difficult  to  use  could  be  potentially  engineered  out  by  redesigning  the 
user  interface.  If  the  input  was  a  graphical  interface,  with  coding  as  a  secondary  method 
of  input,  the  efficiency  of  data  input  could  sky-rocket.  To  draw  the  simulated  eye-model 
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in  a  CAD  program  would  take  a  few  hours  and,  if  cells,  tally  points  and  volumes 
mirrored  a  graphical  rendering,  it  would  eliminate  the  need  to  scrutinize  many 
sections  of  code  after  making  one  minor  change.  While  the  MCNP  Visual  Editor  is 
an  attempt  at  a  graphical  interface,  and  has  proven  itself  a  wonderful  tool  for  verifying 
geometry,  there  is  still  much  to  be  desired  in  the  program.  The  inflexibility  of  the 
current  MCNPX  code  might  easily  translate  into  lost  time  in  dose  planning  and  overall 
treatment  time  in  a  medical  facility.  The  usefulness  of  the  current  laborious  interface 
might  be  questionable  for  day-to-day  practice  in  a  treatment  facility. 

On  a  positive  note,  the  MCNPX  code  performed  surprisingly  well  once  the  program 
was  working  correctly.  The  simulation  run  times  ranged  from  10-13  minutes.  In  that 
time  range  the  MCNPX  statistical  indices  yielded  outstanding  results.  The  relative  error 
was  extremely  low  for  the  typical  and  worst  case  scenarios — 0.008  and  0.0012, 
respectively.  Likewise,  all  the  others  statistical  indicators  resulted  in  values  far  beyond 
the  threshold  for  credible  results.  This  provided  assurance  that  once  the  geometry,  cells, 
and  tallies  were  input,  errors  in  the  code  were  identified  and  corrected,  one  can  expect 
very  precise  results. 

The  MCNPX  code  also  has  additional  helpful  features  that  could  be  incorporated  into 
future  proton  beam  simulations.  For  example,  MCNPX  Version  2.5.0  features  mesh 
tally  plots  which  apply  to  fluxes,  heating,  doses  or  other  tally  quantities  plotted  on  a  3-D 
mesh  surface  (rectangular,  cylindrical,  or  spherical)  independent  of  the  problem 
geometry.  This  could  be  used  to  graphically  illustrate  tissue  heating  effects  due  to 
proton  treatment,  or  doses  over  the  volumes  in,  or  adjacent  to  the  eye.  The  propagation 
of  heat,  dose,  or  other  tally  values  might  be  a  fruitful  study  for  tumors  near  sensitive 
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organs  such  as  the  optic  nerve — while  flux  values  taken  over  the  eye  might  more  clearly 
define  scattering  effects  and  adjacent  dose  distribution  during  treatment. 

Eye  Model  Conclusions 

Greater  detail  could  be  incoiporated  into  the  current  model  of  the  eye — effectively 
expanding  the  types  of  cancerous  tumors  which  might  be  modeled.  Regions  outside  the 
eye  were  neglected  in  the  creation  of  this  model.  Adjacent  organ  dose  is  inferred  by 
adjusting  the  penetration  depth  of  the  beam  with  a  range  modifier,  and  confirming  the 
dose  by  Bragg  peak  behavior  in  the  tallies.  This  may  not  be  an  adequate  in  all  cases. 
The  boney  structure  of  the  skull  for  example,  could  be  a  limiting  factor  in  treatment 
angles  of  the  proton  beam — if  the  treatment  angle  is  too  great,  the  bone  protruding  from 
below  or  above  the  eye  socket  might  interfere  with  treatment.  Likewise,  the  soft  tissues 
surrounding  the  eye  will  receive  dose  when  the  beam  is  placed  outside  transpupillary 
treatment  angles.  Adding  these  structures  would  lead  to  better  understanding  of  doses 
outside  the  eye,  and  potential  limitations  of  treatment  angles. 

Inside  the  eye,  greater  accuracy  could  also  be  attained  by  adding  more  detail.  For 
example  the  cilliary  bodies  and  the  iris  surrounding  the  lens  are  areas  where  cancerous 
growth  can  occur.  With  the  current  eye  model  it  was  not  be  possible  to  simulate  these 
types  of  tumors  and  their  treatments.  Likewise,  greater  resolution  could  be  obtained  by 
differentiating  between  organs  located  closely  together  in  the  eye,  like  the  sclera, 
choroid,  and  the  retina.  For  purposes  of  the  typical  treatment  and  worst  case  scenarios, 
combining  them  into  a  general  outer  structure  of  the  eye  was  sufficient,  but 
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differentiating  might  be  helpful  for  studying  treatment  scenarios  for  a  range  of  tumors  at 
various  stages  of  progression. 
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APPENDIX  B 

Dosimetric  Cells  Volume  Calculation 

The  MCNPX  code  has  limited  ability  to  compute  the  volumes  of  defined  cells. 

The  dosimetric  volumes  in  the  outer  wall  of  the  eye  are  an  example  of  geometry  too 
complex  for  MCNPX.  Thus,  the  volumes  had  to  be  computed  by  hand,  and  input 
directly  into  the  program.  Equation  1  was  used  for  calculating  the  volumes  of 
overlapping  spheres. 

ft(R  +  r  - d)2 (d2  +2dr -3r2  +  2dR  +  6rR-3R2) 

Y2d  ( 

where  R  and  r  are  the  radii  of  the  larger  and  smaller  spheres,  respectively,  and  d  is  the 
distance  between  the  centers.  The  overlapping  volumes  were  calculated  starting  with  the 
section  closest  to  the  optic  nerve,  then  working  successively  forward.  The  previous 
volumes  and  the  aqueous  humor  portion  of  the  eye  were  subtracted  from  the  volume  of 
the  new  dosimetric  cell.  The  results  are  shown  in  Table  4. 


Table  4.  Dosimetric  Cell  Volumes 


Dosimetric  Cells 


Volume  (cm 


29  and  39 
28  and  38 
27and  37 
26  and  36 
25  and  35 
24  and  34 
23  and  33 
22  and  32 
21  and  31 
20  and  30 


0.0090 

0.0322 

0.0760 

0.0789 

0.0822 

0.0863 

0.0914 

0.0978 

0.1059 

0.1166 


MCNPX  Typical  Case  Output  File 
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this  tally  meets  the  statistical  criteria  used  to  form  confidence  intervals:  check  the  tally  fluctuation  chart  to  verify, 
the  results  in  other  bins  associated  with  this  tally  may  not  meet  these  statistical  criteria. 


interval (1,2, 3  sigma):  4.4964E+00  to  4.S077E+00;  4.4908E+00  to  4.S133E+00;  4.4852E+00  to  4.S189E+00 
interval (1,2,3  sigma):  4.4964E+00  to  4.5077E+00;  4.4908E+00  to  4.5133E+00;  4.4852E+00  to  4.5189E+00 
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